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Abstract 



CO , A nonlinear stochastic Schrodinger equation for pure states describing non-Markovian 

' diffusion of quantum trajectories and compatible with non-Markovian master equations is 

presented. This provides an unravelling of the evolution of any quantum system coupled 
^ ' to a finite or infinite number of harmonic oscillators without any approximation. Its power 

fS| , is illustrated by several examples, including measurement-like situations, dissipation, and 

' quantum Brownian motion. Some examples treat this environment phenomenologically as 

an infinite reservoir with fluctuations of arbitrary correlation. In other examples the envi- 
ronment consists of a finite number of oscillators. In such a quasi-periodic case we show 
, the reversible decay of a 'Schrodinger cat' state. Finally, our description of open systems is 

^\ ' compatible with different positions of the 'Heisenberg cut' between system and environment. 



1 Introduction 

> ' 
C , 

, In quantum mechanics, a mixed state, represented by a density matrix pt, describes both an 
ensemble of pure states and the (reduced) state of a system entangled with some other system, 
here consistently called 'the environment'. In both cases the time evolution of pt is given by a 
linear map 

^ : Pt = ^tpo, (1) 



which describes the generally non-Markovian evolution of the system under consideration. Such 
equations describe both an open system in interaction with infinite reservoirs with an arbitrary 
correlation, or a system entangled with a finite environment. In almost all cases, the general 
eq. (H) cannot be solved analytically. Even numerical simulation is most often beyond today's 
algorithms and computer capacities, and thus, the solution of (0) remains a challenge. 

In the Markov limit, eq. (||) simplifies and reduces to a master equation of Lindblad form 

^pt = -i[H, Pt] + ^J2 {i^rnPt, Ll] + [Lra, PiL'^J) , (2) 

where H is the system's Hamiltonian and the operators Lm describe the effect of the environment 
in the Markov approximation. This approximation is often very useful because it is valid for 
many physically relevant situations and because analytical or numerical solutions can be found. 

In recent years, a breakthrough in solving the Markovian master equation (|2|) has been made 
through the discovery of stochastic unravellings of the density operator dynamics. An unravelling 
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is a stochastic Schrodinger equation for states \ipt{z)), driven by a certain noise Zf such that the 
mean of the solutions of the stochastic equation equals the density operator 



Pt = M[\Mz))mz)\\. (3) 

Here M[. . .] denotes the ensemble mean value over the classical noise zt according to a certain 
distribution functional P{z). 

The simplest stochastic Schrodinger equations unravelling the density matrix evolution are 
linear and do not preserve the norm of ipt{z). Such an unravelling is merely a mathematical rela- 
tion. To be truly useful, one should derive unravellings in terms of the corresponding normalized 
states 

\\Mz)\\ 

where now relation (^) can be interpreted as an unravelling of the mixed state pt into an ensemble 
of pure states. Of course, using the normalized states iptiz) requires a change of the distribution 
P{z) — > Pt{z) in order to ensure the correct ensemble mean, with 

Pt{z)^\\Mz)fP{z) (5) 

so that the eq.(^) remains valid for the normalized solutions, 

pt = Mt[\Mz)){Mz)\]. (6) 

We refer to this change (|5|) of the probability measure as a Girsanov transformation |Q] - 
other authors refer to "cooking the probability" or to "raw and physical ensembles" Q, or to 
" a priori and a posteriori states" |Q] . 

In the case of Markovian master equations of Lindblad form (|2|), several such unravellings 
(|6|) are known. Some unravellings involve jumps at random times, others have continuous solu- 
tions. The Monte-Carlo wave function method |^], sometimes called quantum jump trajectories 
1^ 1^, is the best known example of the first class, whereas the Quantum State Diffusion (QSD) 
unravelling is typical of the second class. All these unravellings have been used extensively 
over recent years, as they provide useful insight into the dynamics of continuously monitored 
(individual) quantum processes |9|, In addition, they provide an efficient tool for the nu- 
merical solution of master equations. It is thus desirable to extend the powerful concept of 
stochastic unravellings to the more general case of non-Markovian evolution. First attempts 



towards this goal using linear equations can be found in |11], other authors have tackled this 



problem by adding fictions modes to the system in such a way as to make the enlarged, hy- 



pothetical system's dynamics Markovian again |12, O, 14]. In our approach, by contrast, the 



system remains as small as possible and thus the corresponding stochastic Schrodinger equation 
becomes genuinely non-Markovian. 

Throughout this paper we assume a normalized initial state ipo{z) = ipQ of the system, 
independent of the noise at t = 0. Such a choice corresponds to a pure initial state po = |V'o)('0o| 
for the quantum ensemble and correspondingly, to a factorized initial state ptot = Po ® Penv of 
the total density operator of system and environment. 

In this article we present for the first time the nonlinear non-Markovian stochastic Schrodinger 
equation that unravels the dynamics of a system interacting with an arbitrary 'environment' of 
harmonic oscillators, finite or infinite in number. For a brief overview of the underlying micro- 
scopic model see Appendix C. In the Markov limit, this unravelling reduces to QSD and will 
therefore be referred to as non-Markovian quantum state diffusion. Our results are based on 
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the linear theories presented in , where the problem of non-Markovian unravellings was 

tackled from two quite different approaches. The linear version of the non-Markovian stochastic 
Schrodinger equation relevant for this paper, unifying these first attempts, was presented in [^] 
for unnormalized states. 

Here we present examples of the corresponding normalized and thus more relevant theory. 
We include cases where the environment is treated phenomenologically, represented by an expo- 
nentially decaying bath correlation function, and cases where the 'environment' consists of only 
a finite, small number of oscillators - in section ^ of even just a single oscillator. The latter case 
corresponds to periodic (or quasi periodic) systems, that is to extreme non-Markovian situations. 
Before presenting examples in sections ^, |^, and ^, all the basic equations are summarized in the 
next section ^ Several open problems are discussed in section 0, while the conclusive section |8| 
summarizes the main achievements. 

2 Basic equations 

In this section we summarize all the basic equations. Let us start by recalling the case of 
Markovian QSD, providing an unravelling of the Markovian Lindblad master equation (P). 

2.1 Markov case 

The linear QSD equation for unnormalized states reads 

^ipt = -iH'ilJt + LiltfO zt-^L'^Lipt, (7) 
where zt is a white complex-valued Wiener process of zero mean and correlations 

M[z;zs] = 6{t-s), M[ztzs]=0, (8) 

and o denotes the Stratonovich product [p!^ ]. 

The solutions of eq. (^) unravel the density matrix evolution according to the master eq. (^) 
through the general relation (^. Here, equation (^) is written for a single Lindblad operator L, 
but it can be straightforwardly generalized by including a sum over all Lindblad operators 
each with an independent complex Wiener process z^. 

The simple linear equation has two drawbacks. First, its physical interpretation is un- 
clear because unnormalized state vectors do not represent pure states. Next, its relevance for 
numerical simulation is severely reduced by the fact that the norm of the solutions 

tends to with probability 1 (and to infinity with probability 0, so that the mean square norm 
is constant). Hence, in practically all numerical simulations of <^ the norm tends to 0, while 
the contribution to the density matrix in ^ is dominated by very rare realizations of the noise 
z. 

Introducing the normalized states (Q) removes both these drawbacks. As a consequence, the 
linear eq.(0) is transformed into a nonlinear equation for iptiz)- In this Markov case, the result 
of Girsanov transforming the noise according to (|5|) and normalizing the state is well known 
||2|, m , it is the following QSD evolution equation for the normalized states 

j^i^t = -iH^t + {L- {L)St o {zt + (Lt)i) - i(LtL - {L^L)t)A, (9) 

where (L) = {^tlLli^t)- This equation is the standard QSD equation for the Markov case written 
as a Stratonovich stochastic equation. Notice that it appears in its Ito version in references |^]. 
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The effect of the Girsanov transformation is the appearance of the shifted noise 



zt + {L')u (10) 

entering (^), where zt is the original process of (^). The effect of the normahzation is the 
subtraction of the operator's expectation values. 

2.2 Non-Markovian case 

In the non-Markovian case, the linear stochastic Schrodinger equation generalizing ^ was de- 



rived in reference ||17|] , it reads 

= -iHijjt + Liptzt - [ a{t, s)^ds. (11) 
at Jo ozg 



It unravels the reduced dynamics of a system coupled to an arbitrary 'environment' of harmonic 



oscillators - see Appendix C for a brief overview. Thus, eq.(ll) represents an unravelling of a 
certain (standard) class of general non-Markovian reduced dynamics as in (||). The structure 
of eq.(p!l|) is very similar to the Markovian linear equation (|^: the isolated system dynamics is 
Schrodinger's equation with some Hamiltonian H. The stochastic influence of the environment 
is described by a complex Gaussian process zt driving the system through the Lindblad operator 
L. While this is a white noise process in the Markov case, here it is a colored process with zero 
mean and correlations 

M[z;zs]=ait,s), M[ztzs]=0, (12) 

where the Hermitian a{t,s) = a*{s,t) is the environment correlation function. Its microscopic 
expression can be found in Appendix C. In this paper, we sometimes but not always adopt a 
phenomenological point of view and will often choose a{t, s) to be an exponential (^ exp(— 7|t — 
s\ — iO(t — s)), decaying on a finite environmental 'memory' time scale 7~^, and oscillating with 
some environmental central frequency Q. The Markov case emerges in the limit 7 — > oo. In the 
most extreme non-Markovian case, when the 'environment' consists of just a single oscillator of 
frequency 0, we have the periodic a{t, s) = ex'p{—iQ{t — s)). Finally, the last term of eq.(0) is 
the non-Markovian generalization of the last term of the Markovian linear QSD eq. (0). This 
term is highly non-trivial and reflects the origin of the difficulties of non-Markovian unravellings. 
One can motivate equation ([Tl|) on several grounds. First, it was originally derived from a 



microscopic system-environment model [|17|. In the original derivation the correlation function 
a{t, s) describes the correlations of environment oscillators with positive frequencies. However, 
as can be seen in Appendix C, any positive definite a{t, s) can formally be obtained from some 
suitably chosen environment that possibly includes negative frequency oscillators (Hamiltonian 
not bounded from below). 



Next, as a second motivation, we sketch a direct proof that eq. (11) defines an evolution 



equation (|T]) for density operators. This ensures that the stochastic equation is compatible 
with the standard description of mixed quantum states |]l^, Let po = J2j Pjli^o^) ii^o^ \ 

be any decomposition of the density operator at the initial time (recall that at time zero 
the system and environment are assumed uncorrelated) . What needs to be proven is that pt 
is a function of po only, where pt = This guarantees that pt does not 

depend on the decomposition of po into a mixture of pure states {iV'o''^)}- For this purpose we 



notice that the solution ■0^ of (11) is analytic in z and is thus independent of z* . Hence we find 



^j^{ipt\ = ^iJ^I^^^!^. Accordingly, the evolution equation of |V't)(V't| is linear: it depends linearly 
on |V'o)(V'o|- Since the mean M is also a linear operation, pt depends linearly on pQ. Finally, the 
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positivity of pt is guaranteed by the existence of a pure state decomposition and its normalization 
follows from the fact that eq. (^) preserves the norm in the mean, M[||'(/;t|p] = const as shown 
in Appendix B. 

Thirdly, another set of motivations for eq. (11) is provided by the numerous examples of 



the next sections of this article and by the fact that, by full analogy with the Markov case, 
there exists a corresponding non-linear equation for normalized states, as will be shown in the 
remainder of this section. 



To summarize, eq.(ll) is the basic equation for non-Markovian linear QSD. The functional 
derivative under the integral indicates that the evolution of the state ij^t at time t is influenced 
by its dependence on the noise Zg at earlier times s. Admittedly, this functional derivative is 
the cause for the difficulty of finding solutions of eq ( [ll|) in the general case, even numerical 
solutions. 

We tackle this problem by noting that the linear equation ( p^ ) may be simplified with the 
Ansatz 

^"^^ d{t,s,z)i;t, (13) 



dz 



s 



where the time and noise dependence of the operator 0{t,s,z) can be determined from the 
consistency condition 

= — V'* (14) 
dt 5zs dzs 



with the linear equation ([ll|). The Ansatz ( |T3[ ) is completely general and hence, once the operator 
0{t, s, z) is known, the linear non-Markovian QSD equation ( pT] ) takes the more appealing form 

—i/jt = -iHtpt + Liptzt - a{t, s)0{t, s, z)dsipf (15) 

at Jo 

We are going to show in the subsequent sections how to determine 0{t,s,z) for many interesting 
and physically relevant examples. In most of these cases, in fact, the operator O turns out to 
be independent of the noise z and takes a simple form. 

Being the non-Markovian generalization, eq.(pT|) or equivalently eq.(^) suffers from the same 
drawbacks as its Markov limit (0): the norm of its solutions tend to with probability 1. And 
the cure will be similar. One introduces the normalized states (^) and substitutes the linear 



stochastic Schrodinger equation (15) by the corresponding non-linear one. Its explicit form can 
be rather involved as will be demonstrated in the following sections. 

The derivation of the desired evolution equation of the normalized states ipt requires two 
steps: Taking into account the Girsanov transformation of the noise (^) and normalization. In 
Appendix B we prove that the non-Markovian Girsanov transformation for the noise probability 
distribution Pt{z) (^) corresponds to a time dependent shift of the original process zt, 



zt 



zt+ / a{t,sy{L^)sds. (16) 







This shift and the normalization of the state ipt results, as shown in Appendix B, in the 
non-linear, non-Markovian QSD equation for the normalized state vectors ipt, which takes the 
ultimate form 

j^i^t = -iHiJt + {L-{L)t)i't~zt (17) 
<t, s) ((Lt - (Lt),)0(t, s, z) - ((Lt - (Lt)i)0(t, s, z))t) ds i^t, 
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where Zt is the shifted noise (|T6[). 

Equation (|l^) is the central result of this paper, the non-Markovian, normalized stochastic 
Schrodinger equation that unravels the reduced dynamics of a system in interaction with an 
arbitrary 'environment' of harmonic oscillators - encoded by the properties of the environment 
correlation function a(t,s). In the following sections we will give many interesting examples of 
this non-Markovian Quantum State Diffusion equation (p!7[). 



3 Spin ^ examples 

In this section we use spin ^ examples to illustrate general methods to solve the non-Markovian 
QSD equations (^) (or (p^)) and (17) respectively. These are generally numerical, though 
sometimes analytical, solutions which illustrate certain new features of non-Markovian QSD, 
unknown in the Markov theory. Throughout this section a denote the Pauli matrices. 



3.1 Measurement-like interaction 

This is the simplest example, hence we present it in some detail. Let H = ^<Jz, L = Xaz 
with A a real number parameterizing the strength of the interaction. The harmonic oscillator 
environment is encoded by its correlation function a{t, s) which is left arbitrary in this section. 
First, in order to eliminate the functional derivative in (o), we assume as an Ansatz 



^ = Aa.Vt, (18) 
dzs 

i.e. we choose 0{t,s,z) = Xaz independent of t,s and z in (^). It is straightforward to show 
that, indeed, this Ansatz is compatible with (|l^, i.e. it solves the fundamental linear equation 
(0). 

The corresponding non-linear, non-Markovian QSD eq.([l7|) for the normalized state ■0* reads 

d ~ .u 
di 



tpt = -i'^azipt + A(cr2 - {az)t)A [zt + \ a{t, s)* {az)sds + ^ ot{t, s)ds (cr^)t^ . (19) 



This equation is the generalization of the Markov QSD equation (^ for general environment 
correlations a{t, s). Notice that, indeed, (|^) reduces to the corresponding Markov QSD equation 
(|9|) in the limit of a delta-correlated environment (one has Jq a{t, s)f{s)ds —>■ ^f{t) for any 
function f{t)). 

Eq. ( |T^ ) shows the effect of the non-Markovian Girsanov transformation (^). It induces not 
only the shifted noise ([l6|), but also leads to an additional shift due to the implicit z^-dependence 
of t(^n as explained in detail in Appendix B. Numerical simulations of (p!9|) are shown below. 

In order to find the reduced density matrix of this model, we solve analytically the linear 
non-Markovian QSD equation (|l5|). Using ([Tsl) we find 



d uj 

-j^tpt = -i-^CTzTpt + X(7ziptzt + a{t, sfdsipf (20) 



From the explicit solution of this equation we obtain the expression for the ensemble mean 

Pii(O) Pi2(0)e-^W 

P2l(0)e-^W* /922(0) 



Pit) = M UMM = ( .^.^^-kty '"l^.n^ ) > (21) 
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with F{t) = iujt + 2)? /q ds Jq du{a{s, u) + a*{s, u)). Taking the time derivative, one can show 
that this density matrix is the solution of the following non-Markovian master equations 

Pt = -iT;[(^z,Pt]- — {a{t,s) + a{t,s))ds [az,[az,pt\] (22) 
2 z Jo 

= -'i'T:Wz,Pt]+ }C{t,s)psds, (23) 

2 Jo 

where the "memory super-operator" IC{t,s) acts as follows on any operator A: 

AC(t,s)yl = -y(a(t,s) + a(t,s)*)e-'^'/»''"r^''W"'^)+-("''')*)[a^ (24) 

Let us now turn to actual simulations of this example. In Figure la we show non-Markovian 
QSD trajectories from solving (|l^) numerically with = 2u! and an exponentially decaying 
environment correlation function a{t,s) = ^exp(— 7|t — s\) with 7 = 0; (solid lines). For this 
exponentially decaying environment correlation function the asymptotical solution is either the 
up state or the down state ((cr^) = ±1), while the ensemble mean M[{az)] remains constant 
(dashed line). Thus, as in the standard Markov QSD case, the two outcomes 'up' or 'down' 
appear with the expected quantum probability: Prob(j;'^Vt = I T)) =I(T IV'o)^- Notice that 
for these non-Markovian situations, the quantum trajectories are far smoother than their white- 
noise counterparts of Markov QSD Q . We emphasize that if the environment consists of only a 
finite number of oscillators, represented by a quasi-periodic correlation function a{t, s), no such 
reduction to an eigenstate will occur. 

In Figure lb we compare the average over 10000 trajectories of the non-Markovian QSD 
equation ( |T7|) with the analytical ensemble mean (|2l|) and see very good agreement. This con- 
firms that indeed, both the memory integrals in eq.(p^) arising from the Girsanov transformation 
of the noise are needed to ensure the correct ensemble mean. 



3.2 Dissipative interaction 

This is the simplest example with a non-self adjoint Lindblad operator. Again we set H = ^ctz, 
but now we choose L = Xa- = X^{ax — i(Jy) describing spin relaxation. Also, the environmental 
correlation function a(t, s) and thus the quantum harmonic oscillator environment can be chosen 
arbitrary. 

First we have to replace the functional derivative in ([TT|), and we try an Ansatz ([l3| ) of the 
form 

^-^ = f{t,s)a.^t, (25) 



with f{t,s) a function to be determined. The consistency condition (|T^ of our Ansatz (25) 
leads to the condition on f{t, s): 

dtf{t,s)a-i^t = [-i^a,-XF{t)a+a.,f{t,s)a-]^l^t (26) 

= {iio + XF{t))f{t,s)a^tl;t (27) 



with 



F{t) = f a{t,s)f{t,s)ds. (28) 
Jo 
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Hence, if a^Tpt / 0, the function /(t, s) must satisfy the following eq. 



dtf{t,s) = {iLO + XF{t))f{t,s) (29) 

with initial condition f{s,s) = X. The corresponding non-Markovian QSD equation (|l^ for 
normalized state vectors V't reads 

i't = -i'^CTzi't - XF{t){a+a^ - {a+a-)t)i't (30) 

+ A((j_ - {a^)t)^Jt (^zt + X a{t, s)* {a+)sds + {a+)tF{t)^ 



with F{t) determined from (28) and ([29|). For a given a{t, s), the non-Markovian QSD equation 
( |30| ) can be solved numerically, having F{t) determined numerically from (p9[). Note that in the 
Markov limit, the correlation function a{t, s) tends to the Dirac function 5{t — s). Consequently, 
F{t) tends to the constant ^f{t,t) = ^ and one recovers the standard Markov QSD eq. (^). 

It turns out that non-Markovian QSD can exhibit remarkable properties, unknown in the 
Markov theory. In order to highlight these features, we proceed analytically and assume expo- 
nentially decaying environment correlations a{t,s) = Then we see from ( |2^ ) 



and (g8|) that the relevant function F{t) in (13 Q ) satisfies 



F{t) = -7F(t) + i{uj - n)F{t) + AF(t)2 + ^ (31) 
with initial condition F{0) = 0. With 7 = 7 — i{uj — Q) the solution reads 



F(t) = i -^^i^Unh (|v¥^+ «ctanh (^^,^)) ^ (32) 

For the remainder of this section we assume exact resonance: = lo and thus 7 = 7. Let us 
first consider the case of short memory or weak coupling, 7 > 2A^. For long times, F(t) tends 

to ^7 — — 27A^^ /(2A). For large 7 this asymptotic value tends to ^, which corresponds to 
the Markov limit (|^), as it should. 

More interesting, let us consider the opposite case of a long memory or strong coupling, 
7 < 2A2. In this case, F{t) diverges to infinity when the time t approaches the critical time 
tc = (^7^ + 2 arctan(7/-\/2A27 — 7^)^ /\/2A^7 — 7^. What happens is that at time tc, the first 
component of the vector Tpt vanishes, hence CF-iptc = and eq.(p9|) no longer holds. Indeed, 
the second term of eq.(|30|) becomes dominant and drives the spin to the ground state in a 
finite time, which we prove below in terms of the density matrix. In Figure 2a (for individual 
trajectories) and Figure 2b (for the ensemble average over 10000 runs) we see this effect from 
solving the non-Markovian QSD equation numerically, where we choose A^ = = a;, so that 
Lvtc = Ivr ~ 4.71. For t > tc the state tpt is constant. This is an example where a stationary 
solution is reached after a finite time! It is the first example of a diffusive stochastic Schrodinger 
equation which is at the same time compatible with the no-signaling constraint (ie the evolution 
of mixed states depends only on the density matrix, not on a particular decomposition into a 
mixture of pure states) and has no "tails" (does not take an infinite time to reach a definite 



state), see the discussions in |21, |2^. In [^] it is proven that such a feature is impossible for 
Markov situations. Notice that this peculiar feature holds at resonance only. 

Finally, we note that for the intermediate case 7 = 2A , one has F{t) = — > A for t ^ 
00, again approaching a constant value (the reader may find it helpful to adopt our convenient 
convention for the choice of units: [zt] = [A] = [f{t)] = [F{t)] = [l/Vt] and [a(t, s)] = [1/t]). 
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In order to determine the corresponding master equation for the reduced density opera- 
tor, we solve the hnear QSD eq. ( [Tl| ) where we make use of the change of variable: (j)t = 
^-i^a^t+Xa+a- F{s)ds^^^ After some computation and taking the ensemble mean analytically, 
one gets 

V P2i(0)e Jo ^1 i-p^^{t) J 

This proves that whenever Re|^ Jq F{s)ds^ diverges for a finite time, the density matrix pt reaches 
the ground state in that finite time and thus all pure state samples have to do so as well. For 
the time evolution of this reduced density matrix one gets 

Pt = -i'^[a,,pt] + XiF{t)+F{ty)(^a^pta+-^{a+a-,pt}^ (34) 

= -i-[crz,Pt] + / JC{t,s)psds, (35) 
^ JO 

where the "memory super-operator" IC{t,s) acts as follows on any operator A: 

/C(t,s)A = -y(a(t,s)+a(t,s)*)( 2e~^/.' ^(")'^V_^a+ - A} (36) 

- 2(e-^ ^(")'^" - l)a+a_ Ac7+a_) 

In Figures 2a and 2b we illustrate this example (A^ = 7 = = w) for exponentially decaying 
correlations. All individual non-Markovian quantum trajectories reach the ground state in the 
critical time tote ~ 4.71 (Fig 2a). Taking the ensemble mean over 10000 trajectories, we find 
very good agreement with the analytical result of the reduced density matrix (Fig. 2b). 



4 More examples 

4.1 Model of energy measurement 



This case, H = L = LK is a straightforward generalization of section |3^. Again, the environment 
correlation a{t, s) can be chosen arbitrary. We find O = H in (13) and the non-Markovian QSD 
equation ( p^ for the normalized states reads 

^Jt = -iH^t - {H^ - {H^)t)i't f a{t, s)ds (37) 
at Jo 

+ {H - {H)t)^t[zt + j\{t,sf{H)sds + j\{t,s)ds {H))j. 
For the corresponding master equation we find 

Pt = - t a{t, s)ds H[H, pt] - f a{t, sYds [pt, H]H, (38) 
Jo Jo 

hence Ti{Hpt) is constant, contrary to the individual expectation values (H)^^. 

The eigenvectors of H are stationary solutions of the non-Markovian QSD eq. (|37[) . Thus, if 
the noise is large enough, all initial states tend asymptotically to such an eigenstate, as in Markov 
QSD. However, if the noise has long memory, as for example in the extreme case of periodic 
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systems (see section such a reduction property clearly does not hold. The exact conditions 
under which eq. ( |37D describes reduction (localization) to eigenstates are not known. Notice 
however that if the correlation decays smoothly such that Jq a{t, s)ds tends for large times t to 
a real constant, and if {H)t converges for large times to a fixed value, then the non-Markovian 
equation (^7|) tends to 

^4>t = -iHi^t - {H^ - {H^)Mtconst + {H - {H)St {zt + const{H)t) . (39) 

The long time solutions of this equation are the same as the long time solutions of the corre- 
sponding Markov approximation. The latter is the Markov QSD equation, hence the asymptotic 
solutions tend to eigenstates of H. The previous subsection |3.1| provides an example of this more 
general statement for H = ^Cz- 



4.2 A simple toy model 



In this subsection we use a simple toy model |24] to illustrate that the non-Markovian QSD 



equation ( p^ may contain unexpected additional terms that cancel in the Markov limit. Con- 
sider H = p and L = q and an arbitrary environment correlation function a{t,s). Then the 
Ansatz ( p!^) fo r replacing the functional derivative with some operator satisfying the consistency 
condition ([l^ ) reads 

^ = (g-(t-.))V^,. (40) 



Thus, the non-Markovian QSD eq. ( p?^ ) takes the form 

—^t = -ipi^t - [q^ - {q^)t) a{t, s)dsi^t (41) 

+ {q- {<l)t)i^t (^zt + a{t, s)*{q)sds + a{t, s)ds{q)t^ 

+ {q-{q)t)'4't[ {t - s)a{t,s)ds. 
Jo 

The first two lines of this non-Markovian QSD equation could have been expected, since they 
have the same form as in the previous examples, see for instance eq. (^). The last line of 
the above equation, however, has no counterpart in the previous examples. Clearly, it vanishes 
in the Markov limit {a{t,s) — > 5{t — s)), when the non-Markovian QSD equation (^) for this 
model reduces to the Markov QSD equation (P). 



4.3 Quantum Brownian motion model 

In this subsection we consider the important case of quantum Brownian motion of a harmonic 



oscillator [25|, that is we choose H = ^{p +q ), L = Xq, and arbitrary environmental correlation 



a{t, s). As shown in Appendix C, the basic linear non-Markovian QSD equation for this quantum 
Brownian motion case is again the fundamental linear equation dll 



It turns out that the functional derivative in (11) is more complicated in this case, because 
0(t,s,z) depends explicitly on the noise z. However, fortunately, this dependence is relatively 
simple. Indeed, let 



— = 0{t,S,Z)'ll)t 

oz. 



f{t,s)q + g{t,s)p + i [ ds'j{t,s,s') 
Jo 



(42) 
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The consistency condition ( [L4D leads to the following equations for the unknown functions 
fit, s), git, s) and jit, s,s') in (|4|): 

dtfit,s) = vogit,s)+ifit,s) t ds'[ait,s')git,s')] (43) 







2igit,s) ds'[ait,s')fit,s')]-i ds'[ait, s')jit, s' , s)] 



dtgit,s) = -ojfit,s)-igit,s) ds'[ait,s')git,s')] (44) 

Jo 

jit,s,t) = git,s) (45) 

dtjit,s,s') = -igit,s) f ds"[ait,s")jit,s",s')]. (46) 



These equations have to be solved together with the non-Markovian QSD equation (|T^. 

If, for simplicity, we assume exponentially decaying environment correlations ait, s) = 
lg-7l*-«l a,nd introducing capital letters for the integrals, Xit) = Jq ait, s)xit, s)ds, for x = 
f,g,j, one obtains the simpler closed set of equations 

Fit) = Y~ + '^'^^^^ ~ - ^^-^(t) (47) 

Git) = -7G(t) - wF(t) - iAG(t)2 (48) 

Jit) = ^Git)-2^Jit)-iXGit)Jit), (49) 

where J(t) = ait, s')Jit, s')ds' . The initial conditions read F(0) = G(0) = J(0) = 0. Finally, 
Jit, s) can be determined from the solutions of the above equations, we get 

Jit,s) = AG(s)e-/^^+'^^(^'))''^'. (50) 
Hence, the non-Markovian QSD equation for quantum Brownian motion becomes 

j^i^t = -iHiPt - iq^ - {q\)^tFit) - iqp - {qp)t - p{q)t + {p)t{q)StGit) (51) 



+ iq - {q)t)4^t { zt + ait,s)*{q)sds + {q)tFit)-i J(t, s')(^s' + / ais' , s)* (q) sds)ds' 
\ Jo Jo Jo J 

Let us make some comments about this non-Markovian QSD equation. First, recall that it 
corresponds to the exact solution of the quantum Brownian motion problem |25] of a harmonic 
oscillator. Next, this example shows a new feature that we didn't encounter in the previous 
examples: the noise Zt enters the equation non-locally in time. Thirdly, terms involving the 
operator qp appear, although there are no such terms either in the Hamiltonian or in the en- 
vironment operator L = Xq. Finally, since this equation is exact, it is a good starting point to 
tackle the quantum Brownian motion problem using this new approach and to find its proper 
Markov limit. In connection with this last point, we emphasize that the master equation cor- 
responding to eq. (|5T|) necessarily preserves positivity |Q because it provides a decomposition 
of the density operator into pure states at all times. However, these questions and numerical 
simulations are left for future work. 
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4.4 Harmonic oscillator at finite temperature 



As another important example of an open quantum system we briefly sketch the case of a 
harmonic oscillator H = uja^a coupled to a finite temperature environment through L_ = A_a. 
As explained in detail in Appendix C, the finite temperature also induces absorption from the 
bath, which has to be described by a second environment operator = A+a^. Hence, the linear 
non-Markovian QSD equation (^) has to be modified and involves two independent noises 
and z. 



t 5 



—ipt = -iHtpt + X-aiptZ-t 



Sip 

6zs 



a+{t,s)P^ds, (52) 
ozT 



see eq. (|9^ ) in Appendix C. This equation can be solved with the following Ansatze: 



5ih 
Szs 

5zt 



f-{t,s)a+ f ds'j-{t,s,s')z+ 
Jo 



/+(t,s)a^+ / ds'j+{t,s,s')z^, 



Ipt- 



(53) 



(54) 



Using similar techniques as in the previous subsection, the evolution equations for f±{t,s) and 
j±{t, s, s') can be obtained and thus the resulting non-Markovian QSD equation can be written 
in closed form. A new feature of this example, again unknown in the Markov case, is that each 
of the two environment operators L_ and is coupled to both noises. 



5 Harmonic oscillator coupled to a few oscillators: decay and 
revival of Schrodinger cat states 

The case of a harmonic oscillator coupled to a finite or infinite number of harmonic oscillators 
all of which are initially in their ground state (zero temperature), H = uja^a, L = Xa, is very 
similar to the damped spin ^ example treated in subsection The Ansatz = f{t,s)aijjt 
similar to ( |25| ) holds with f{t,s) and F{t) satisfying the same equations (p9| ) and (p8|). Thus, 
the non-Markovian QSD equation (^) for this situation reads 

^Jt = -iuja^aiPt+{a-{a)t)Mzt+ a*{t,s){a^)sds+XF{t){a^)t)-XF{t){a^a-{a^a)t)ipf (55) 
dt Jo 

Again, this non-Markovian QSD equation reduces to the Markov equation (^) for a{t, s) = 



6{t — s) since in this case F{t) = 4 according to (28). As in the case of a dissipative spin 



(subsection for exponentially decaying bath correlations at resonance, the system oscillator 
may reach its ground state in a finite time, provided the correlation time I/7 is long enough. 

Notice also that (^) preserves coherent states \f3). The time evolution of the complex number 

labeling these coherent states is given by 

$t = {-iio-F{t))f3t. (56) 

More interesting than a coherent state initial condition is the case of a superposition — /?) 
of two symmetric coherent states, known as a 'Schrodinger cat' |^]. If the correlation decays, 
so does the 'cat'. If, in contrast, the environment consists only of a finite number of oscillators, 
then the 'cat' will first decay, due to the localization property of QSD, but since the entire 
system is quasi-periodic, the 'cat' will then revive! 
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As an illustration, we simulate the extreme case where the 'environment' consists of only 
a single oscillator. It thus models the decay and revival of a field cat state in a cavity that is 
isolated from the outside, but coupled to a second cavity, to which it may decay reversibly. Such 
an experiment on reversible decoherence was proposed recently in |28]. In this simple case, the 
environment correlation function reads 

a{t,s) = e~'^^'-'\ (57) 

where is the frequency of the single 'environment' oscillator. Figure 3 shows the time evolution 
of the Q-function of such a 'Schrodinger cat' in phase space for Q = 0.5w and a coupling strength 
between the two oscillators of O.loj. Apart from an overall oscillatory motion due to the 'system' 
Hamiltonian ioa^a, we see how the cat first decays but later becomes alive again. Further 
investigations of stochastic state vector descriptions of such reversible decoherence processes are 
left for future investigations. It is worth mentioning that depending on the stochastic process, 
the cat my subsequently decay into either of its two components. 

6 Shifting the system- environment boundary 

In this section we consider a situation where the 'Heisenberg cut' between the system and the 
environment is not obvious. Since the non-Markovian QSD equation provides the exact solution 
of the total system-environment dynamics, the description of the system does not depend on 
this cut. This is in contrast to the usual Markov approximation, where the position of the cut 
is crucial. As an example, let us consider a system consisting of one spin ^ and one harmonic 
oscillator, the two subsystems being linearly coupled. Assume moreover that the spin ^ is 
coupled to a heat bath at zero temperature, see figure 4. The total Hamiltonian reads: 

Hiatal = Hi + H2 + H12 + Henv + Hj (58) 

with 



Hi 


= 0", 

2 


(59) 


H2 
H12 

Henv 


= x{(^-a^ + o-+a) 


(60) 
(61) 
(62) 


Hi 




(63) 



We can either consider the spin-oscillator system coupled to a heat bath, or consider only 
the spin coupled to a heat bath and coupled to an auxiliary oscillator, as illustrated in Figure 
4. In the first case, we can consider the Markov QSD description, ie a family of spin-oscillator 
state vectors il^t{C) indexed by the complex Wiener processes ^j. In the second case, using non- 
Markovian QSD we have a family of spin ^ state vectors (t>t{S,,z) indexed by the same plus 
the non-Markovian noise zt with correlations 

M[z-*Zs]=e-'''^^^'-'\ (64) 



The (linear) stochastic eqs.(ll) governing ipt and (pt read 



ijjt = -i{Hi + H2 + Hi2)iJt + X(^-ipS - —(J+a-i^t (65) 
= -iHi^t + XcT-Mt - ^(T+^-<Pt + X(^-(l^tZt - t e-''"2(*-^)^ds (66) 

Z Jo OZs 
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where A is a function of the Xa;'s, that is of the strength of the spin-heat-bath coupling. 

A natural question in the present framework is to study the " Heisenberg cut" : compare the 
states of the spin ^ averaged over the noise z with the mixed state obtained by tracing out the 
second oscillator (Tr2) from the 1-oscillator-spin states, i.e. we ask whether the equality 

M,[|c/>i(e,z))(<Ai(e,z)|] =rr2(|Vt(6)(^t(e)l) (67) 

holds. According to the general non-Markovian QSD theory presented in this paper, the spin^ 
state should be independent of the position of the Heisenberg cut. Below we illustrate this 
feature using the present example. 

By assumption the oscillator starts in the ground state: V'o = (/"o ® |0). Hence, the state '^t 
can be expanded as 

V't = co(t)| i,0)+Ci(t)| T,0)+C2(t)| 

where 

Co = A^(t)ci ^-y Co 

ci = -(^Y + y)ci - iXC2 (70) 
C2 = -i((a;2 - y)c2 + xci)- (71) 



Tracing out the single harmonic oscillator, one obtains the spin ^ state (in the jj, basis) 

PI = Tr^{\Umm)\) = f K , isl'u ,2) ■ (72) 

V CqCi |Co| + |C2| / 



We now turn to the alternative description of the same situation, but with the 'cut' between 
the spin ^ and the oscillator. In order to solve eq. (66) we make the usual Ansatz 

^ = /(t,s)cT_,/.t, (73) 

where the consistency condition ( |l^ ) leads to dtf{t, s) = {iuJi + ^+xF{t))f{t, s), where f{t, t) = 
X and F{t) = Jq a{t, s)f{t, s)ds. Consequently, 

F{t) =x + {iooi -iuJ2 + Y + xF{t))F{t). (74) 
Using the notations (pt = ^^o(OI i) + ^i(*)| T) one gets 

^'0 = « Y^'o + (A6 + XZt)vi (75) 

VI = -(iy + y + Xi^(i))t'l. (76) 

Note that since vi is independent of zt, vi{t) is itself independent of z, hence, 

^M, [vo] = i yM^ [vo] + X^tvi . (77) 



Averaging over the z-noise, one obtains the spin ^ state (in the jj, basis) 

« - M. [I0,(f, = {J^. . (78) 
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Finally, a straightforward comparison of eqs. (6S,70,71) and (74, 75,^) shows that cq = -/^^[fo], 
ci = vi and C2 = —iFvi. Hence, 3 of the 4 entries of the matrices pi and p2 are equal. The 
equality of the fourth entry follows from the general feature that linear non-Mar kovian QSD 
preserves the mean of the square norm. 

This completes the proof that pi = P2- the spin ^ state is independent of the position of the 
Heisenberg cut, for all times and all realizations of the heat bath induced noise This illustrates 
the general fact that non-Markovian QSD attributes stochastic pure states to systems in a way 
which depends on the position of the Heisenberg cut, but which is consistent for all possible 
choices of the cut. See Fig 4. for the illustration of these relationships. This is in opposition to 
the case prevailing in Mar kovian unravellings. 

7 Open problems 

This paper is the first presentation of non-Markovian QSD. Admittedly, there remain many 
open questions and a lot of work has still to be done to exploit all the possibilities opened up 
by this new approach. In this section we list some of the open problems: 

1. The ultimate goal would be to develop a general purpose numerical simulation program. 
However, at present no general recipe is known. 

2. When do the long time limit and the Markov limit commute? A question which is of 
particular interest for quantum Brownian motion. 

3. If the initial condition is not factorized, the present approach must be generalized. 

4. In the Markov case unravellings exist both with continuous trajectories and with quantum 
jumps. In the non-Markovian case, the only unravelling known at present is the continuous 
non-Markovian QSD described in this article. What about non-Markovian unravellings 
with quantum jumps? 

5. In the Markov case, continuous QSD unravellings exist for real or pure imaginary noise, 
as well as for complex noise. What about the non-Markovian case? It seems that in the 
present case complex noise is essential. 

6. Note that most of the non-Markovian master equations used in this article have known 
analytical solutions. In these cases, the general Zwanzig form [pg| ] of the master equation: 

pt = f lC{t-s)psds (79) 





with the memory kernel IC{t — s) could be rewritten as a Lindblad type master equation 
with time-dependent coefficients. Then, the master equation can also be simulated using 
Markov QSD with time dependent coefficients. However, if the solution of the master 
equation is not known explicitly, or does not lead to a Lindblad type equation, then 
numerical simulation has to use the non-Markovian QSD theory. It would be interesting 
to illustrate non-Markovian QSD for more of such examples and to study the conditions 
under which a non-Markovian problem can be treated with Markovian unravellings. 

How does non-Markovian QSD compare with consistent histories [^] and other approaches? 
For instance, it was shown in |31] that the solutions of the non-Markovian eq. ( [T7| ) can be 



considered as conditional states in the framework of a " hybrid" representation of the fully 
quantized microscopic system. 
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8. What is the perturbation expansion of the non-Markovian QSD eq. (p^) in terms of the 
memory time 7~^? The zeroth order term would be the Markov QSD eq. (^), what about 
the higher orders? 

9. Finahy, non-Markovian QSD should be applied to open problems in physics, where non- 



Markovian effects are relevant, such as semi-conductor lasers [12|, or atom lasers [p^ ]. 

8 Conclusion 

We present a stochastic equation for pure states describing non-Markovian quantum state dif- 
fusion, compatible with non-Markovian master equations. We illustrate its power with several 
examples. In essence, we show that quantum (finite or infinite) harmonic oscillator environments 
can be modeled by classical, complex Gaussian processes, entering the non-linear, non-Markovian 
stochastic Schrodinger equation for the 'system' state which we derive in this paper. 

Several authors have proposed stochastic pure-state descriptions of such non-Markovian sit- 
uations using fictions modes added to the system in such a way as to make to dynamics of the 
enlarged hypothetical system Markovian |12, 13]. Others, [O] treat a non-Markovian problem 



with an explicitly time-dependent Markov unravelling. In our approach, by contrast, there are 
no additional modes, hence the system is as small as possible, and the stochastic Schrodinger 
equation becomes genuinely non-Markovian. This is of interest for efficient numerical simula- 
tion and high-focus insight into the relevant physical processes. Also, non-Markovian quantum 
trajectories are in general much smoother than those of Markov processes, which might even 
help to reduce further the numerical effort. 

Let us stress an important conceptual difference between Markov QSD and non-Markovian 
QSD. In the Markov case, one starts from a master equation for mixed states and associates to 
it a stochastic Schrodinger equation. The master equation may either be derived from a micro- 
scopic model, or merely be based on phenomenological motivations In the non-Markovian 
case, on the contrary, one starts from the stochastic Schrodinger eq. (|TT1). The existence of a 
master equation is guaranteed by the microscopic model summarized in Appendix C. In general, 
however, the explicit form of this master equation is not known. Nevertheless, this existence 



ensures that the corresponding stochastic Schrodinger equation for normalized states (17) does 
not allow arbitrary fast signaling, despite its nonlinearity |20| ]. 

From a pragmatic point of view, the Hamiltonian and environment operators in eq.([l^) can 
either be derived from a microscopic theory, or be merely based on phenomenological motiva- 
tions. Non-Markovian master equations are almost always exceedingly difficult to treat, even 
numerically. However, one can always start from the non-Markovian QSD approach of this 
paper, which appears thus more fundamental than the master equation approach. 
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A Frequency representation 

It is sometimes useful to express the noise by frequency components z^j'. 

zt = Y.z^e^\ (80) 

where the frequencies a; can take positive as well as negative values. Also the correlation function 
can be written in Fourier representation: 

a(t,s) = a(t - s) = ^a<^e-''^(*-'), > 0. (81) 

The correlation of the Fourier components of the noise is trivial: M[2;*za] = ^ld\oi-w In this 
representation the distribution functional becomes a simple Gaussian distribution over all z^'s: 

P[z) = M ^^J-YP^ (82) 

and the states V't become functions of the frequency amplitudes z^^ of the noise. We can then 
write the fundamental linear non-Markovian QSD eq. (|Tl|) in terms of them: 



^J,t = -iHiJt + 5](Le*"*z^ - Lta,e-*"*/-)Vt. (83) 

at cij oz^^ ^ 

This frequency representation is a helpful tool to discuss the mathematical properties of the non- 
Markovian stochastic Schrodinger equation (|l^), as we do in Appendices B and C. Remember 
that in eq. (|83| ) we assume the initial condition to be independent of the noise: iPq{z) = i/jq- 



B Girsanov transformation for non-Markovian QSD 



As time goes by, Girsanov transformation distorts the distribution P{z) ( p^ ) of the complex 
noise into Pt{z) according to eq.(^). In frequency representation, we have 



\z 12 



P,(z)=mMz)fe^p(-Y.^\. (84) 

We assume that at t = the state ipo is normalized and does not depend on z. So, initially, 
Po{z) is identical with P{z). 

We find the time evolution of Pt{z) from the linear non-Markovian Schrodinger equation ([l] 
in frequency representation (^3|). Using eq.(|8^), we find 



^^P,iz)=M{Mz)\j^M^))exp (-E (85) 



Now we make a crucial observation. The solution ipt{z) of (p3|), with initial condition ipt{z) = ipo, 
is analytic in all z^j's. Then it follows that dliptiz)) /dz* = d{'4)t{z)\/ dz^^ = for all z^. Hence, 
when inserting eq. ( |8^ ) into eq. (|85[), we can substitute 



{^,{z)\L^ ^Uz)) = ^{L^)tmz)f, (86) 

OZi.) OZi,j 
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and we obtain 



dt 



-Ptiz) 



E 



-iujt 



d_ 

dz,,^ 



{L^)tPtiz) + c.c. 



(87) 



This is a remarkable result. It shows that the Girsanov transformation is equivalent to a drift 
of the random variable z. We read off the drift velocities directly from eq. 



d_ 

dt 



-iujt 



One can see that the Girsanov transformation preserves the normalization of the distribution 
Pt{z). This has the immediate consequence that the non-Mar kovian stochastic Schrodinger 
equation (11) preserves the mean norm of the quantum state: 



P{z)dz = / Pt{z)dz = 1. 



M 



Now we are going to derive the stochastic non-Markovian Schrodinger equation for the nor- 
malized states tptiz) = il^tiz)/\\iptiz)\\j where ipti^) is the unnormalized solution of the linear 
stochastic equation ([ll|). First, we solve the drift eq. ( |88| ) for the trajectories Za,{t), with the 
initial conditions -2(^(0) = z^ for all uj: 



Zuj + 



(90) 



where {L'^)t = {'4't{z{t))\L'^\4't{z{t))) / {4't{z{t))\'ipt{z{t))) . The Girsanov-transformation (^) leaves 
invariant the probability of the noise z along the above trajectories: 



Pt{z{t))dz{t) = P{z)dz 
for all z^. Hence, we can write the stochastic unravelling (P) as follows: 



pt = Mt \Mz))mz)\ 



M 



\Mm)){Mm)\ 



(91) 



(92) 



The mean value on the very right refers to the simple undistorted distribution P{z). To calculate 
it, one has to express %l)t{z{t)) as a function of the initial amplitudes z^ = Z(^(0). Remember that 
^t{z) is the solution of the linear non-Markovian equation (11) or (^) with initial condition 
il>t{z) = ipo- The additional time dependence of ipt{z{t)) through z{t) appends a new term to 
the evolution equation of these 'Girsanov-shifted' states, so that we find the following stochastic 
evolution equation: 



Ql^t + l^Z^—^t 



dZn 



iHi^t + V Le'"*z^ - (Lt - (Lt)i) /* a{t, s)d{t, s, z)dsA 



(93) 



(94) 



where we used (|^), ([l5|), and (^8|). Finally, these states have to be normalized. The resulting 
evolution equation for the normalized states is our central result, given by eq. (p^ . In the 
time domain, the shifted noise (pO|) takes the form (16). 
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C Review of the linear non-Markovian theory 



Here we briefly review the microscopic origin of tlie linear non-Markovian stochastic Schrodinger 
equation - see [15, 16, 17 1. The linear non-Markovian QSD equation results from a standard 



model of a system interacting with an environment of harmonic oscillators, represented by a set 
of bosonic annihilation and creation operators a^^,alj. The interaction term Hj between system 
and environment is chosen to be linear in the a^jS and arbitrary in the system operator L: 
Hj = Xuj{Lalj + L^a^), with some coupling constants Xuj- Thus, the model is defined by 

Htot = Hgys + Hi + Henv (95) 
= Hsys + ^Xu;{Lal + L^Uuj) +^ujala^ (96) 

Solving this total closed system in a clever way leads to the linear non-Markovian stochastic 
Schrodinger equation ( [ll| ) for the system state il^t{z). As initial condition we assume a factorized 
form ptot = |V'o)(V'o| ®Pt for the total density operator, with all bath oscillators initially in some 
thermal state pT = ®ujPui{T). 



C.l Zero temperature 

In it was shown that if all the environment oscillators are initially in their ground state 
(T = 0), the dynamics of the reduced density operator pt = tr envPtot{t) of the model (|95|) can 
be unraveled {pt = M[\ipt{z)) {'iptiz)\]) using the linear stochastic Schrodinger equation (pl|), 

^V't = -iHx/jt + Li^tzt - f a{t, s)^ds (97) 
at Jo ozg 

where the colored complex stochastic processes zt with zero mean satisfy 

M[z*tz,] =Y.xl e~^"(*^^) ^ a(t, s), M[ztZs\ = 0. (98) 

UJ 

We see the microscopic origin of the bath correlation function a(t, s) at zero temperature. For 
real physical systems we have a; > in (|98[). To model an arbitrary time-translation invariant 
correlation function, one needs environment oscillators with negative frequencies as well. 



C.2 Finite temperature 

In order to derive the linear non-Markovian QSD equation at finite temperatures, we use a 
simple mathematical trick, well known in field theory |3^: the non-zero temperature density 
operator pT of the heat bath can be canonically mapped onto the zero-temperature density 
operator (the vacuum) of a larger (hypothetical) environment. The problem at T > is thus 
reduced to the problem at T = 0, whose linear non-Markovian QSD equation (|97|) we already 
know. The resulting finite temperature linear non-Markovian QSD equation is 

^i^t = -iHi/jt + LtPtzi - [\-{t,s)^ds + L^i;tzt - L f a+{t,s)^ds. (99) 
at Jo ozs Jo ozj 

It thus depends on two independent processes z^^zf with zero means and with temperature 
dependent correlations 

M[zt*z-] = J2{n^ + l)xl (*-'^) ee a-{t,s), M[z^ z'] = (100) 
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and 

Mlz^zt] =Y.^u.xl e^"(*-^) = a+{t,s), M[z+z+] = 0. (101) 

Here, n^j = (exp ^ — 1)~^ denotes the average thermal nmiiber of quanta in the mode uj. We 
identify these terms as describing the stimulated (n) and spontaneous (+1) emissions {Lz~) and 
the stimulated absorptions (n) from the bath (L^^z^). Notice also that for T ^ 0, all the n^^ 
tend to zero and (p9|) reduces to (|97|), as it should. 



C.3 Finite temperature and L = 

In the case of a selfadjoint coupling operator L = = K, the finite temperature result can 
be simplified considerably by introducing the sum process zt = -\- z^' having zero mean and 
correlations 

M[z1zs] = a'^{t, s) + a^{t, s) = a{t, s) = ^ xtiC^^uj + 1) cosa;(t — s) — ismuj{t — s)] 

UJ 

M[ztzs] = 0. (102) 

Notice that (2na; + l) = coth so that a{t, s) is nothing but the well-known bath correlation 

kernel of the so-called quantum Brownian motion model [^] . In terms of this single process Zt , 
the linear non-Markovian QSD equation at finite temperature (^9|) takes the simple form of the 
zero-temperature equation (|97[) involving just one noise zt 



-Jt = -iH^Pt + K^jtzt -K I a{t, s)^ds, (103) 
at Jo ozs 

with the temperature dependent a{t,s) of (102). For K = q the position operator, this un- 



ravelling was first introduced in [|16| , derived from the exact Feynman- Vernon path integral 
propagator of this model. 
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Figure captions 



Fig la: Quantum trajectories of the non-Markovian QSD equation for the 'measurement '-hke 
case H = ^az, L = Xoz and an exponentially decaying bath correlation function a(t,,s) = 
^ exp(— 7|t — s|). We choose = 2a;, 7 = 0; and an initial state |V'o) = (1 + 2i)| t) + (1 + 01 i)- 
Displayed is the expectation value (cjz) of several solutions of the non-Markovian QSD equation 
(solid lines) and their ensemble average (dashed line). 

Fig lb: Same parameters as in Fig la: Here we compare the ensemble average of the Bloch 
vector using 10000 quantum trajectories of non-Markovian QSD (solid lines), with the analytical 
result (dashed lines). 

Fig 2a: Quantum trajectories of the non-Markovian QSD equation for the dissipative case 
-f^ = ^Cz, i = Aa_ and an exponentially decaying bath correlation function a(t, s) = ^ exp(— 7|t— 
s\ — iQ{t — s)). We choose = 7 = cj and resonance Q = to. As initial state we use 
IV'o) = 3| T) + I !)• Displayed is the expectation value (cr^) of several solutions of the non- 
Markovian QSD equation (solid lines) and their ensemble average (dashed line). At the finite 
time Lotc = Itt 4.71, all individual trajectories reach the ground state. 

Fig 2b: Same parameters as in Fig 2a: Here we compare the ensemble average of the Bloch 
vector using 10000 quantum trajectories of non-Markovian QSD (solid lines), with the analytical 
result (dashed lines). 

Fig 3: Reversible decay of an initial Schrodinger cat state \ipo) = \a) + \ — a) with a = 2. 
The contour plots show the Q-function of a non-Markovian quantum trajectory of a harmonic 
oscillator (u), coupled to just a single 'environment' oscillator ($7 = 0.5a;), initially in its ground 
state. The coupling strength between the two oscillators is 0.1a;, and the time step between two 
successive plots is 2.27/a;. 

Fig 4: Shifting the 'system-environment' boundary. First, we consider the 'spin - single 

oscillator' system with state iptiCi^ coupled to a heat bath with noise Alternatively, we can 
consider the 'spin' only as the 'system' ^^(^,2;), coupled to the 'single oscillator -|- heat bath' 
environment (noises ^tj^t)- In non-Markovian QSD, both descriptions are possible and lead to 
the same reduced spin state. 
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